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We present the exact solution of two Heisenberg-liquid models of particles with arbitrary spin 
S interacting via a hyperbolic long-range potential. In one model the spin-spin coupling has the 
simple antiferromagnetic Heisenberg exchange form, while for the other model the interaction is of 
the ferromagnetic Babujian-Takhatajan type. It is found that the Bethe ansatz equations of these 
models have a similar structure to that of the Babujian-Takhatajan spin chain. We also conjecture 
the integrability of a third new spin-lattice model with long-range interaction. 



I. INTRODUCTION 

Spin fluctuations are a common feature of strongly correlated electron systems such as high Tc superconductors 
and heavy fermion compounds. Exchange interactions of the Heisenberg type, originally introduced to describe the 
properties of insulating magnets, have now been realized playing a central role in doped Mott insulators and heavy 
fermion systems. Two related models are the t — J model and the Kondo lattice model. In these systems, the spin 
exchange among the itinerant electrons or between the conduction electrons and the local moments induce fascinating 
physical phenomena, such as very large mass enhancements and non- Fermi liquid behavior. For some low-carrier- 
density systems, especially in low dimensions, the small number of carriers and their reduced mobility do not provide 
an effective mechanism for screening and long-range interactions should be considered. A special class of Integra 
systems with long-range interactions are the Calogero-Sutherland model (CSM)ld'cl and its SU{N) aeneralizationsJ 
The lattice versions of these models, which generically aije referred to as Haldane-Shastry models,|j can be obtained 
by freezing the orbital dynamics within a proper scheme.l3 In this letter, we present two integrable models of particles 
with arbitrary spin interacting via a Calogero-Sutherland potential with a hyperbolic space dependence. In case I the 
spin-spin coupling is of the simple Heisenberg exchange form, i.e. Si ■ Sj, while in case II, the interaction is of the 
Babujian-Takhatajan model (BTM) type.B Finally, based on the structure of the Bethe ansatz equations (BAE), -an 
S'[/(2)-invariant high-spin chain model, which corresponds to the long-range-coupling generalization of the BTM,eI is 
conjectured to be integrable. 

The structure of the present paper is the following: In the subsequent section, we construct the,modcl Hamiltonians. 
The energy spectra and the BAE are derived via the so-called asymptotic Bethe ansatz (ABA).a In sect. HI, wc study 
the thermodynamics as well as the ground state properties. Concluding remarks are given in sect. IV. 

II. MODELS AND BETHE ANSATZ 

The general form of the model Hamiltonian we shall consider is 

where Xj is the position of the particle j, 7 > is a real constant, and gij is a spin-dependent coupling, which for our 
two cases is: 

=2S,-S, + 2SiS+l), (2) 

for model / and 

5y = ~4SQ2s{S^ ■ S,) = [2S - hj] [2S+1~ f„], (3) 

for model //, where 
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which has eigenvalues 0, • • • , 25 and Q2S is just the Babujian polynomial^ Hence, in case / the spin coupling is a 
simple antiferromagnetic exchange, while in case // it is ferromagnetic and of the Babujian- Takhatajan form. It is 
interesting to compare model / to the t — J model. Their physics is quite similar, although in the present model the 
coordinate Xj is continuous while it is restricted to a lattice in the t — J model, and that here the interaction is of 
variable range. For some high Tc compounds, the next-nearest-neighbor spin exchange is believed to be relevant. In 
the present model the range of the interaction is governed by the free parameter . For instance, for sufficiently 
large 7 the interaction is short-ranged, while the long-range Calogero-Sutherland interaction is readily recovered in 
the limit 7^0. We note that in both cases, I and II, the total spin is conserved in the two-body scattering processes 
and the Hamiltonian (|l]) has a basic S'Li/2)-invariance. 

The asymptotic Bethe ansatz (ABA)H is a powerful tool to derive the energy spectrum for models with nonlocal 
interactions, like e.g. the CSM. If the model is integrable it is sufficient to know the asymptotic long-distance behavior, 
i.e. the scattering phase shifts, which can be obtained without the full knowledge of the wave functions. The ABA 
also assumes that the particles are asymptotically free, i.e., that they do not form charge bound states. The ABA only 
needs the two-body phase shifts (or two-body scattering matrix), but requires an independent proof of the integrability 
of the model Hamiltonian. The integrability of the present model is related to the unitary matrix model of Minahan 
and Polychronakos,u who have shown that the CSM with a coupling strength = JijJji is exactly soluble. In Ref. [ 
^ Jij was chosen to be the SU{q) angular momentum. For the 2S'-fold SU{2) case, we can map JijJji to (St + Sj)"^, 
which is of more interest to condensed matter problems and corresponds to our case /. In fact, a potential J^/r^ 
can be generated naturally from the dynamics of two free particles in three dimensions in spherical coordinates. Such 
a potential describes the angular dynamics of the center of mass of two free particles. If the angular momentum is 
considered a "spin" , the problem is reduced to a one-dimensional one and the scattering matrix must be factorizable, 
which already infers the integrability of many particles in one dimension with such a potential. Hence, for case / 
our choice of gij is qwite natural. On the other hand, case // can be realized by a proper quantization of Jij using 
projection operators.cl 

To derive the two-body S'-matrix, we note that the two-particle scattering occurs in 25 + 1 different channels of 
total spin I ranging from to 25. In a given channel of spin /, the coupling strength is l{l + 1) for case / and 
(25 — 1){2S — I + 1) for case //. For convenience we will express our results in terms of the operator kj defined above. 
Below we study the two cases separately. 



A. Model / 

The two-body transmission 5-matrix can be derived following Sutherland's method.i Consider the two-particle 
case. The asymptotic wave function for jxi — X2I — > 00 can be written as 

*(a;i,a;2) ^ Ai2(12)e*'^^^i+*'^^^^ Ai2(21)e''=^^i+''=i^% (4) 

for xi < X2 and 

^{xi,X2) A2i{2l)e'^^='^+'^^''^ + A2i{l2)e'^^''^+'^^''\ (5) 

for xi > X2, where fci, ^2 are the asymptotic momenta carried by the particles and the A' s are constant coefficients, 
which determine the 5-matrices. For example, 5i^2(fci — ^2) = ^2i(21)/Ai2(12). By solving the Schrodinger equation 
for two particles we obtain 



q (h h\- r(l + z%)r(l-H/,,- -^fc,,-) 



where kij — {ki — kj)/2^, r(a;) is the gamma function and Pij is the spin permutation operator. In fact, Sij is 

a polynomial of Si ■ of order 25 and is proportional to the Lax-operator of the BTM.eI Hence, it satisfies the 
Yang-Baxter relationE2l SijSi^mSj.m — Sj,mSi,mSij . For example, when 5 — 1/2, we have 



_ % +.(1/2 + 25. ■5j-) 
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while for S' = 1 it is 



^ 2Si ■ Sj{l — ikij + Si ■ Sj) — kfj + ikij 2 

~ (z% + l)(i%+2) • ^ ' 

Because the system is non-diffractive, the wave function must asymptotically be given by plane waves ^{x\Q) ^ 
J2p^{-P\Q) ^^P['''J2^=i^Qj^Pj] when \xi — Xj\ — > oo for all pairs where xq-^ < • • ■ < xqj^, and P and Q are 
permutations of (1, 2, • • • , iV). With the periodic-boundary conditions of the wave function, the momenta kj are 
determined by a set of Yang eigenvalue equationsEj 

• • • %i%Ar . . . %,+iVl/(P|Q) = e-^^^'^^iPlQ). (9) 

For convenience, we introduce an auxiliary spin and define Sqj — <S'oj(A — kj), so that SojiO) — Poj. Further 
we define the (25* + 1) x (25* + 1) monodromy matrix 

2o(A) = So,j ■ ■ ■ S'cA'S'o,! • ■ ■ S'oj+i. (10) 

By summing over the auxiliary spin of the matrix To(A), it is easy to show that the trace 

ti"oTo{kj) — Sj,j-i ■ ■ ■ Sj,iSj,]\i ■ ■ ■ Sjj+i 

is just the transfer matrix of the BTA/I^ with an extra multiplying factor and lattice disorder (we consider a contimuuB. 
model). The eigenvalue problem Eq.(9) can be solved following the standard method of the algebraic Bcthc ansatzptil 



e 



»fc,L - J If TT fcj - fc» - 2m7 I yr kj - Xg + 2iSj 
J-J- J--'- kj-h + 2m7 [ kj - Xa- 2iSj' 



n—l i.i^j 



-p- Aa - kj - 2iS^ _ yr ^g-^p- '2il , 
A„ + 2^57 "^;J,J„A„-A;3 + 2*7' ^ ^ 

where Xa are the rapidities of the spin waves and M is the total number of spin wave quanta. Here we dropped 
a multiplying factor exp[«7r(2S' + 1)A^] in the first of Eqs.(ll), because it is irrelevant in the thermodynamic limit. 
The first set of Eqs. (11) is determined by the eigenvalue of the transfer matrix troTo^kj) — exp(— ifcjL), while the 
second set arises as the condition that ensures that the Bethe states are eigenstates of the Hamiltonian. The energy 
eigenvalues are given by 

N 

E = J2kl (12) 
i=i 

by acting the Hamiltonian on the asymptotic wave function. 



B. Model II 



For our second model the two-body S'-matrix is given by 



^ _ r(i + z%)r(i + 25-f,,--i%) ^ 



With the same procedure, we obtain the BAE as 

M 

n 



e 



ikiL TT ^3 -^a 2iSj 



kj — Xa + 2iSj 



-pr Xg - kj - 2iS"f ^ yr Xg - Xp - 2ij 

f}^Xg-k,+2iS^ ^li^A„-A^ + 2z7- ^ ^ 



The energy eigenvalues are still given by Eq.(12). 
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III. THERMODYNAMICS 



A. model I 



Since gij > 0, the potential is repulsive and no charge bound state can exist. This can be seen clearly in the two- 
particle case. In fact, the two-particle wave function can be solved exactly in terms of two hypergeometric functions 
which has no bound state solution in the parameter region we considered. The solutions oLthe BAE, Eqs.(ll), are 
therefore classified into real charge rapidities kj and strings of length n of spin-rapidities.E3 In the thermodynamic 
limit L 00, N oo, N/L Ue, we denote with p{k), (T„(A) and Ph{k), tT^(A) the rapidity densities and the 
densities of their holes, respectively. By taking the thermodynamic limit of the BAE, we have 

^ = p{k) + ph{k) + I3'p{k) - J2 B2S,n<Jn{k), 

n 

oo 

^nW + ^m,n<JniX) = S2S,„p(A), (15) 

m=l 

where the integral operators are given by B' — B2S,n = [25* + n + 1 — 2Z], A^^n — [m + n] + 

2[m -I- n — 2] + • • • + 2[|m — n| -f- 2] -I- [\ni — and [n] is the integral operator with the kernel 

The free energy is expressed as 

F/L = f{k^-p- Sh)p{k)dk + ^ n/i / an{X)d\ 

-T j[{p{k) + ph{k)) ln(p(fc) + phik)) - p{k) Inpik) - ph{k) \nph{k)]dk (17) 
-TY.1 [(^n(A) +a:j(A))ln(a„(A) -f fTi'(A)) - a„(A)lna„(A) - a':,{\)\na':,{X)]d\. 

n •' 

where /i, T and h are the chemical potential, the temperature and the magnetic field, respectively. It is convenient to 
introduce the quantities ({k) = adk) / p{k) and ?7„(A) = cr^(A)/cr„(A). By minimizing the free energy with respect to 
the particle dendities, we obtainllj 

InC = 1 + B' ln(l + C') " E ^2S,„ ln(l + (18) 

ln(l + ?7„) ^ + 51 ^":» li^(l + ^rV/) - ^2S,n ln(l + C"^)- 

m— 1 

Substituting the above equations into Eq.(17), we obtain 

F/L^^—j \n[l + C\k)]dk- (19) 

The thermodynamic BAE (18) can be further simplified. For convenience, we define G — [l]/([0] -I- [2]). With the 
following operator identities 

^l,m GA2.m = Sl.m, 
Bm,n — G[Bm,n+l + Bm.n-l] — <5m,nG, (20) 
Bl^m -— Gi?2,m — Sl^mG, 

we can rewrite the second equation of Eqs.(18) as 
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Inrjn = -(5„,2sGln(l + C"') + G[ln(l + r?„+i) + ln(l + rjn-i)], (21) 

with the boundary conditions 770 = 0, and hm^^oo Inr/n/n = h/T = 2xo. Prom Eqs.(20) obviously GAm,n = -Bm,n, 
so that from the second of Eqs.(18) we have 

Gln(l + 7?2S) = + E ^rn,2S + ) - ^^25,25 ln(l + C"')- (22) 

m 

Substituting the above equation into the first equation of Eqs.(18), we get 

InC = + B' ln(l + C"^) - G[S2S,2sln(l + C"') + ln(l + ms)] 

= ^^+B"Hl + C')-Gln{l+ms), (23) 

where B" = GE^liPm + 1]. 

For T ^ and h wc obtain from Eq.(23) that C ^ for fc < y/JI and C ^ oc' for A; > ^/Jl. This implies that 
the charges form a Fermi sea. In addition, in the same hmit 772s while all other rjn tend to some constant as 
indicated by Eq.(21). Therefore, the spin quanta form a Fermi sea of 2S'-strings in the ground state, as in the case 
of the BTM. The ground state properties are then given in terms of p{k) and <J2sW (which describe the two Fermi 
seas). The rapidities are densely distributed in the intervals [—K,K] and [—A, A], respectively, where in zero-field 
A —> 00 for finite 7, and the integral equations are 

P{k) = r ^'C" - k')p{k')dk' + r B2s,2s{k - A)a25(A)dA, (24) 

J_K J-A 

i-K 



I A2s.2s(A-A')a2s(A')dA'= / B2s,2s{^- k)p{k)dk, 

J-A J-K 



where B'{k), B,„,n(k) and ^,n,n(A) arc the kernels of B' , Bm.n and Am.n, respectively. The elementary excitations 
are uniquely determined by the dressed energies ec{k) (of charges) and es(A) (of spins), which for h ^ satisfy 



B'{k-k')ec{k')dk' + / B25,2s(fc- A)e,(A)rfA, (25) 

-K J-A 

/A 
A2s,2s(A - A')e«(A')rfA' = / B2s,2s(A - k)e,{k)dk, 
-A J-K 

The group velocities of the two bands are Vc = e'^{K)/2'Kp{K) and Vg = e'g(A)/27rcr2s(A), where the prime denotes 
derivative with respect to the rapidity. By integrating the second of Eqs. (25) for /i = (A = 00) wc have = 2Af, 
i.e. as expected the ground state is a spin singlet. At low T the system then is a two-component Luttinger liquid 
with charge dynamics of central charge 1 and spin dynamics with central charge 25. The latter is very similar to the 
BTM. 

In the Calogero-Sutherland limit, all the above quantities can be derived explicitely, since for 7^0, the integral 
operator [m] becomes a (5-function, 



ee(fc)=2e,(fc) = ^^, (26) 
Vc = Vs = 27rne(6' + 1), 

where 9{x) is the step function and the Fermi momentum \s K = wUeiS + 1) with tie = N/L. Without magnetic 
field the particles occupy each spin state equally, so that the distribution function rig (A;) {s = —S, - ■ ■ , S) \s that of an 
exclusion statistics with the statistical weight g = {S + 1){2S + 1), 

_ e{K-\k\) 

"«('^^- (5 + l)(25+l)- 
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B. Model II 



The solutions of the BAE for model // are still given by real kj and A-strings of length n. In the thermodynamic 
limit, the BAE, Eqs.(14), yield 

1 °° 

— = p{k) + ph{k) + Y,Bn,2S<Jn{k), (28) 



ji=i 



im,«Cr„(A) + a';,{X) = i?2S.n/5(A). 

The free energy is still given by Eqs.(17) and (19). Using a similar procedure, i.e., by minimizing the the free energy 
fimctional with respect to p{k) and cr„(A), we readily obtain the thermal BAE of model // 

lnC = ^B2S,nln(l + r?-^), 

h ^ 

ln(l + ry„) = ^ + ^ i^,„ ln(l + r;„i) + Bss.n ln(l + C"'), (29) 



or equivalently 



In C = (fc2 - /i)/r - G ln(l + r^2s) - GB^sas Ml + C"'), (30) 
lnr;„ = (5„,2sGln(l + C"^) + Gln[(l + 77„+i)(l + Vn-i)], 

with the same boundary conditions as in Eq.(21). A strikingly different feature of this model is that when T — > 0, 
/i — > 0, all rjn — > oo, which indicates that there are no flipped spins in the ground state. Hence, the ground state is 
ferromagnetic, since all the electrons align along a spin direction. The rapidity distributions in the ground state are 
then 

p{k)^^e{K-\k\), (31) 
<yn{X) = 0, 

where K is the cutoff (integration limit or Fermi momentum) given by 

K = TTOe. (32) 

The group velocity of the charge excitations at the Fermi surface is 

Vc — 27rne. (33) 

The spin sector is no longer a Luttinger liquid, but has a quantum critical point. A single spin-wave excitation with 
momentum p, has an excitation energy es(p) ~ if p is small. At low but finite T, the charge sector behaves almost 
like non-interacting spinless fermions, while for the spin sector we have a ferromagnetic Babujian-Takhatajan spin 
liquid, so that asymptotically as T — > the specific heat and the spin susceptibility arell3 

G^ri, x-T^^ln^^T. (34) 

The low-T specific heat can also be inferred from the spin-wave dispersion relation. 

In the Calogero-Sutherland limit 7^0, bothm, the charge and spin sectors, obey ideal fractional statistics. The 
thermal BAE can then be solved analytically, since in this limit the integral equations reduce to algebraic ones. 
Solving for r]n separately for n < 2S and n > 2S we obtain 

,„^ ""'[^+^^"] -l, n<2S, 
sm a 

smh'[in-2S)xo + (3] 

Vn = -2 1: n> 2S, (35) 

smh xo 
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where a and /3 are functions of C(fc) determined by 

sin[(2S' + l)Q;]sinha::o = sinh/3sina, 

sin(25a)v/TTC^= ^^|sinh(/3-2;o)|. (36) 
smli a;o 

The solution of Eqs. (35) and (36) yields ??„ as a function of (. Substituting the result into the first of Eqs.(30) we 
obtain a relation for ({k), which uniquely determines the free energy. The distribution function for the charges is 

n(fc) = -|-ln[l + r^(fc)]. (37) 

A careful calculation also shows that the logarithm in Eq.(34) disappears, i.e., x ~ T^^. This is due to the free nature 
of the particles in the limit 7 — s- 0. 



IV. CONCLUDING REMARKS 



In conclusion, we constructed two integrable models with long-range spin couplings. In model /, the coupling 
between spins is of the antiferromagnetic exchange form, while in model //, the spin coupling is ferromagnetic of 
the Babujian-Takhatajan type. The BAE for the spin sector are the same as for the BTM (here the kj introduce a 
lattice disorder). As we have learned from ordinary integrable models, almost all lattice models have a continuum 
counterpart in the sense that their BAE sharej-a common structure. Examples are the spin-1/2 isotropic Heiseaberg 
chainEj with the (5-potential bos©H gas modclJiS the Hubbard modelEZl with the J-potentiaL Fermi gas model,E3 and 
the SU (A^)-invariant t — J modelEJ with the SU (A^)-invariant J-potoitial Fermi gas modelJlj However, no continuum 
counterpart to the BTM with contact interaction was known so far.E3 The present work fills this gap in the integrable 
model family. 



Finally, we conjecture that the following lattice model 

N 



H = sinh2 7 y . " • ^"^ , 
sinh [7(771 — n)] 



(38) 



is integrable for the following reasons: (a) With the proper freezing process (eliminating the chargjs dynamics), the 
continuum models may be reduced to Eq . (^8|) , since the BAE have the same form as for the BTMO (b) For ^ 00, 
Eq.(^8|) flows to the BTM, which has been demonstrated to be integrable. 
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China. 
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